NPS-OR-94-012 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


Approved  for  public  release;  distribution  is  unlimited. 
Prepared  for: 

Ballistic  Missile  Defense  Office 
Washington,  D.C.  20301 

19941201  113 


NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943-5000 


Rear  Admiral  T.  A.  Mercer  Harrison  Shull 

Superintendent  Provost 

This  report  was  prepared  for  and  funded  by  the  Ballistic  Missile  Defense  Office, 
Washington,  DC. 

Reproduction  of  all  or  part  of  this  report  is  authorized. 

This  report  was  prepared  by: 


rr.  KwOX 

KNEALE  T.  MARSHALU^^^ 
Professor  of  Operations  Resbarch 


Reviewed  by;  Released  by: 


Department  of  Operations  Research 


REPORT  DOCUMENTATION  PAGE 

Fom  Approved 

OUB  No.  0704-0188 

Pidtbcteponing  burden  iortMColoaian  of  nbmwlian  is  estimaied  ID  average  1  hour  per  response,  sidudng  tie  time  far  revrawringinstruciions,  searching  exisling  data  souroes, 
gejherjrig  and  mainairiirig  tie  daarieedod,  arid  cotnplelirig  arid  reviewfing  tie  rialotlicin  of  irMomiaton.Seridcernrtiari»regnrdngts8  burden  eelirnalB  or  ariyoiner  aspect  of  fiis 
eeloction  of  iiloiTnaliori.(rictKfcig  suggestions  for  reducing  tiis  burden.  10  Washington  HeadquarlBrsServioe8.Di»BdDralB  for  WbfrrialionOperalioris  arid  Reports,  1218  Jeflersw 
Oavis  Highway,  Suite  1204,  Arwigton,  VA  22202-4302,  and  to  tie  Office  of  Mariagenient  and  Budget,  Paperwork  Reduction  Project  (070^188),  Washinaion,  DC  20S03 

1 .  AGENCY  USE  ONLY  (leev*  bluA)  2.  REPORT  DATE  3.  IMPORT  TYPE  AND  DATES  COVERED 

September  1994  Technical 

4.  TITLE  AND  SUBTITLE 

The  Roles  of  Coimterforce  and  Active  Defense  in  Countering 
Theater  Ballistic  Missiles 

5.  FUNDING  NUMBERS 

6.  AUTHOR(S) 

Kneale  T.  Marshall 

7.  PERFORMING  ORGANIZATION  NAME(S}  AND  ADORESS(ES) 

Naval  Postgraduate  School 

Monterey,  CA  93943 

B.  PERFORMING  ORGAMZATION 

REPORT  NUMBER 

NPS-OR-94-012 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S}  AND  AOORESS(ES) 

Ballistic  Missile  Defense  Office 

Washington,  DC  20301 

10.  SPONSORING  /  MONITORING 

AGENCY  REPORT  NUMBER 

1 1.  SUPPLEMENTARY  NOTES 

1 2a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

12b.  DISTRIBUTION  CODE 

1 3.  ABSTRACT  (Maximum  200  words) 

This  report  contains  the  formulation  and  analysis  of  a  model  to  measure,  compare,  and 
contrast  the  effects  of  counterforce  (pre-launch  attack)  and  active  defense  (post-launch  attack) 
against  tactical  ballistic  missiles  (TBM's).  It  is  shown  that  without  counterforce  an  active 
defense  system  could  require  an  impractical  number  of  weapons  to  counter  incoming  missiles 
and/or  their  warheads.  This  number  is  shown  to  decrease  geometrically  as  effective 
counterforce  is  used,  so  that  the  expected  number  of  warheads  killed  increases  dramatically 
with  counterforce  that  is  only  modestly  effective.  Actual  distributions  of  warheads  reaching  the 
target  area  are  shown  to  be  complex  mixtures  of  binomial  distributions.  It  is  shown  that  normal 
approximations  to  these  distributions  based  on  the  easily-calculated  means  and  variances  often 
agree  poorly  with  the  actual  distributions.  This  is  especially  true  when  using  effective 
counterforce. 

14.  SUBJECT  TERMS 

search,  road  search,  theater  ballistic  missiles  (TBM),  missile  laimchers 

15.  NUMBER  OF  PAGES 

27 

16.  PRICE  CODE 

1 7.  SECURITY  CLASSIFICATION  1 8.  SECURITY  CLASSIFICATION  1 9.  SECURITY  CLASSIFICATION 

OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

Unclassified  Unclassified  Unclassified 

20.  LIMITATION  OF  ABSTRACT 

UL 

I"- .  . . . . .  ''  — '  I  - -  I  ... 

I  NSN  7540>01-280>5S00  Standard  Fonn  298  (Rev.  2-89) 


Prescribed  by  ANSI  Std.  239-18 


The  Roles  of  Counterforce  and  Active  Defense  in 
Countering  Theater  Ballistic  Missiles 


Kneale  T.  Marshall 
Professor  of  Operations  Research 
Naval  PostgraduateSchool 
Monterey,  CA  93943 


September,  1994 


Accesion  For 

,4 _ 

NTIS 

CRA&! 

DTIC 

TAB 

□ 

Unani'.oun.ced 

□ 

Justification 

1 

By  . 

Disc  ibutic;i  1 

1 

A' 

^3;iCjbA‘i 

Ibb!!! . 

Avci!  a 

ndlor 

D’iSt 

Ci2i 

A-/ 

L... 

ABSTRACT 

This  report  contains  the  formulation  and  analysis  of  a  model  to  measme,  compare,  and  contrast 
the  effects  of  counterforce  (pre-launch  attack)  and  active  defense  (post-launch  attack)  against  tac¬ 
tical  ballistic  missiles  (TBM’s).  It  is  shown  that  without  counterforce  an  active  defense  system 
could  require  an  impractical  number  of  weapons  to  counter  incoming  missiles  and/or  their  war¬ 
heads.  This  number  is  shown  to  decrease  geometrically  as  effective  counterforce  is  used,  so  that 
the  expected  number  of  warheads  killed  increases  dramatically  with  counterforce  that  is  only  mod¬ 
estly  effective.  Actual  distributions  of  warheads  reaching  the  target  area  are  shown  to  be  complex 
mixtures  of  binomial  distributions.  It  is  shown  that  normal  approximations  to  these  distributions 
based  on  the  easily-calculated  means  and  variances  often  agree  poorly  with  the  actual  distributions. 
This  is  especially  true  when  using  effective  counterforce. 


1  bitroduction 

TTie  purpose  of  this  paper  is  to  present  and  analyze  a  model  of  theater  ballistic  missile  CTBM) 
launcher  and  missile  flight  operations  so  that  con^arisons  can  be  made  of  the  effectiveness  of  var¬ 
ious  strategies  to  counter  the  threat  The  model  presented  here  extends  earlier  analysis  and  results 
found  in  Conner,  Ehlers,  and  Marshall  [1993]. 

Figure  1  shows  a  schematic  of  the  operations  of  a  TBM  launcher  and  the  missile  assumed  in  this 
report  Launchers  are  expected  to  be  stored  in  some  fixed  storage  area.  When  hostilities  are  about 


Figure  1:  Schematic  of  Theater  Ballistic  Missile  Operations 


to  commence  the  launchers  will  move  to  a  forward  area  for  assembly,  fueling  and  mating  with  the 
missiles.  From  there  a  launcher  will  move  to  its  launch  area,  and  after  launch  will  return  to  the  for¬ 
ward  area  to  prepare  for  the  next  launch.  We  assume  that  each  launcher  has  the  potential  to  launch 
m  missiles,  after  which  it  must  be  taken  out  of  service  for  an  extended  time.  The  reason  could  be 
fliat  it  must  undergo  extensive  repair  and  refit,  or  it  could  run  out  of  missiles.  We  also  assume  that 
each  missile  has  n  (^1)  warheads. 

In  this  papa:  we  assume  that  thoe  are  five  phase  in  the  TBM  oppration  when  the  missile  system 
could  be  attacked.  These  are 
(a)  Counterforce 

1.  Attack  of  the  launcher  with  mated  missile  before  launch  between  assembly  area 
and  launch  site.^ 
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2.  Attack  of  the  launcher  after  missile  launch  either  at  the  launch  site  or  on  return  to 

assembly  area. 

(b)  Active  Defense 

3.  Attack  of  the  missile  during  the  boost  phase, 

4.  Attack  the  missile  on  reentry  before  multiple  wariieads  separate, 

5.  Attack  each  wariiead  in  flie  terminal  phase. 

The  effectiveness  of  attacking  the  system  in  each  of  these  five  phases  is  assume  to  be  summa¬ 
rized  by  a  kill  probability  />,•  for  the  i-th  phase,  or  equivalently  by  a  sxuvival  probabihty  where 
Qi  =  1  -  Pi-  Although  it  is  more  usual  to  formulate  a  model  in  terms  of  kill  probabilities,  survival 
probabilities  are  used  in  Section  2  because  of  the  simplification  that  results  in  ntKxlel  development 
and  presentation  of  results.  Our  objectives  are  to  find  the  probability  distribution,  mean,  and  vari¬ 
ance  of  the  number  of  warheads  reaching  the  target  area  from  each  launcher,  and  the  e7q>ected  num¬ 
ber  of  weapons  required  in  each  phase,  in  terms  of  the  maximum  number  missiles  per  launcher  (wi), 
the  number  of  warheads  per  missile  (n)  and  the  five  survival  probabilities  for  the  five  phases  as 
shown  in  Figure  1.  Using  expressions  for  these  quantities,  in  Section  3  we  connate  the  effect  of 
changing  the  model  parameters  to  demonstrate  that  counterforce,  with  effectiveness  measured  by 
qi  and  will  almost  surely  be  a  necessary  part  of  a  layered  defense  system;  without  at  least  a  mod¬ 
est  success  rate  in  prosecuting  the  launchers  effective  active  defense  may  not  be  feasible. 

2  The  Anti-TBM  Model 

We  build  the  mathematical  model  in  stages  following  the  missile’s  path  from  being  mounted  on 
the  launcher  to  its  or  its  launcher’s  destruction,  or  the  arrival  of  its  warheads  in  the  target  area.  First 
we  develop  the  probability  distribution,  mean  and  variance  of  the  number  of  missiles  that  are  suc¬ 
cessfully  launched  fix)m  a  given  launcher.  These  clearly  will  depend  on  the  counterforce  effort 
against  the  launcher.  Next  we  derive  the  probability  distribution,  mean  and  variance  of  the  number 
of  missiles  that  survive  the  boost  and  reentry  phases.  Finally  we  find  expressions  for  the  probability 
distribution,  mean  and  variance  of  the  number  of  warheads  that  survive  the  final  phase.  The  distri¬ 
bution  of  the  warheads  surviving  to  reach  the  target  area  is  a  complex  mixture  of  Binomial  proba¬ 
bilities.  The  section  ends  with  numerical  examples  to  illustrate  the  results.  A  detailed  analysis  using 
the  model  is  presented  in  Section  3. 

2.1  Launcher  Movement  Phases 

Let  X  be  the  number  of  missiles  launched  from  a  given  launcher  before  it  is  either  destroyed  or 
has  launched  m  missiles.  We  assume  independent  attacks  each  time  the  launcher  attempts  an  out¬ 
ward  journey  to  the  launch  site,  and  similarly  for  each  time  it  attempts  an  return  journey  to  reload. 
Thus  X  is  a  random  variable  that  can  take  on  any  integer  value  from  0  (the  launcher  is  destroyed 
on  the  first  outward  journey)  to  m  (all  attempts  to  destroy  the  launcher  fail).  Note  thatX  >  /  if  and 
only  if  the  launcher  survives  the  first  outward  journey,  and  then  survives  i  succeeding  cycles  back 
to  the  reload  point  and  out  again  to  the  launch  site.  Thus 


^  It  is  understood  that  a  launcher  may  employ  a  numbo'  of  tactics  on  its  way  to  or  fircrni  the  launch  site,  such  as  stop¬ 
ping  in  hide  sites.  The  model  summarizes  the  effects  of  these  strategies  in  a  single  survival  or  kill  {vobability. 
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Pr{X>0}=9i 

Pr{X>l}=«i(^i^2) 

Pr{X>2]  =q,<.giq^ 

•  •  * 

Wxn-1)  =«i(«ii72)"-' 

Pr(X>m)  =  0. 

By  summing  this  cumulative  tail  distribution  (see  result  1  of  the  Appendix): 


m- 1 

£[X]  = 

i<l 


gl(l-(glg2)”) 


(1) 


TMs  equation  holds  if  both  0  ^  <  1  and  0<q2<h  and  is  equal  to  m  when  both  and  ^2  are 

equal  to  1  (zero  effect  in  killing  the  launcher  before  or  after  launch). 

To  find  its  variance  we  need  to  find  its  second  moment  Using  result  2  of  Appendix,  the  second 
moment  of  X  is 


£[X'] 


(m-1)  »i(l-  Wig;)”) 


(2) 


when  both  0  ^  qi  <  1  and  0  ^  ^2  <  1»  and  is  equal  to  when  both  qi  and  qi  are  equal  to  1. 

We  find  the  variance  of  X  in  the  usual  way  by  subtracting  the  square  of  Equation  (1)  from 
Equation  (2). 

We  now  turn  to  finding  the  expected  number  of  weapons  required  in  the  first  two  phases.  Before 
attempting  to  do  this  it  is  necessary  to  make  two  important  assumptions  which  are  assumed  to  hold 
in  all  five  phases.  First,  we  assume  that  every  time  there  is  an  opportunity  to  attack  the  launcher, 
the  missile,  or  one  of  its  warheads,  this  opportunity  is  taken  and  prosecuted  with  a  single  weapon. 
It  may  be  that  in  practise  more  than  one  weapon  is  used,  so  that  the  numb^  determined  by  the 
model  in  this  report  can  be  thought  of  as  lower  bounds.  Second,  the  extreme  case  of  some  kill  prob¬ 
ability  being  zero  in  a  given  phase  can  be  obtained  in  one  of  two  ways,  either  (i)  by  not  attempting 
an  attack  during  that  phase,  or  (ii)  by  attacking  with  a  completely  ineffective  weapon  systenL  Bi 
this  paper  we  assume  that  the  first  of  these  is  true;  any  time  we  use  a  Pi  of  zero  (q,-  of  one)  in  phase 
I  we  assume  no  weapons  are  expended  in  phase  i.  The  expected  numbers  of  weapons  required 
should  not  be  interpreted  as  estimates  of  weapons  requirements  in  actual  operations.  In  this  paper 
they  are  intended  as  an  aid  in  gaining  insight  when  conq>aring  the  effectiveness  of  changing  kill 
probabilities  in  the  various  phases. 

Let  Wbl  and  be  the  numbers  of  weapons  used  in  the  “before  launch”  and  “after  launch” 
phases  respectively  against  the  launcher.  Notice  that  if  no  missiles  are  launched,  zero  (the 
launcher  was  destroyed  on  its  first  outward  journey).  It  is  easy  to  show  that  no  matter  how  many 
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missiles  arc  launched  from  a  given  launcher,  Wj^ = X  and  its  first  two  moments  arc  given  by  Equa¬ 
tions  (1)  and  (2). 

By  following  the  cycle  of  the  launcher  one  can  see  that  the  cumulative  tail  distribution  of 
is  given  by 

Pr{W^L  >  i)  =  if  t  =  0, 2, 4, . . (w-1), 

=  0  if  i  ^  m. 

Using  result  1  of  the  Appendix, 


l-(gig2)” 


and  by  comparing  this  with  Equation  (1)  we  can  see  that 

E[1VbJ=£[X]/<7i.  (3) 

As  our  analysis  progrcsses  through  the  boost  and  reentry  phases,  expressions  arc  found  that  re¬ 
quire  the  probability  mass  function  (pmf)  of  X.  From  the  cumulative  tail  distribution  above  this  is 
seen  to  be 

Px(0)  =  1  -  9i. 


Px(f)  =  ^i(l  ■  ^i^2)(4'i^2y  ‘  4  =  1, 2, . . .,  m-1,  (4) 

Px(m)  =  qi(qiq2T'. 

22  The  Boost  and  Reentry  Phases 

The  boost  phase  and  reentry  phase  survival  probabilities  arc  q^  and  q^  respectively  (see 
Figure  1).  Let  the  number  of  missiles  surviving  both  of  these  phases  (per  launcher)  be  Y.  Clearly 
this  is  also  a  random  variable,  and  if  we  assume  that  the  attempt  to  shoot  down  a  given  missile  in 
either  phase  is  independent  of  the  outcomes  of  earlier  or  later  atten5)ts  at  othw  missiles,  the  con- 
ditioni  random  variable  [ITX]  has  a  Binomial  distribution  with  parameters  X  and  q^q^.  Thus 
£[yiX]  and  Var\Y\X\  -Xq-^q^  (1  -  q^q^.  By  unconditioning  on  X,  the  erq)ected  number  of 

warheads  surviving  the  through  the  reentry  phase  is 

£[n=?3?4£m  (5) 


where  E\X\  is  given  by  Equation  (1). 

The  variance  of  T  is  found  using  the  standard  conditional  variance  argument, 

Vca-m  =  fxEWirrnX]]  +  Varx[E[nX]], 
so 

Var\Y\  =  ^3^4  (1  -  93?4)£IX]  +  {qiq^^Var[X\,  (6) 

where  Equations  (1)  and  (2)  are  used  to  find  l^PQ. 
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To  find  the  pmf  of  Y,  note  that 


Pnx(/10  =  *0;t.93?4) 

where  0^;<i,0<^3^4<l,  and  b(j',i,P)  -  • 

Unconditioning  on  X  we  find 

m 

Pr  0)  =  X*  ^ ^3^4^ Px  (*)  *  i  =  0*  2 . "*• 

i^j 

where  the  Px(0’s  are  given  in  Equation  (4). 

Let  Wb  and  Wr  be  the  number  of  wea5)ons  used  in  the  boost  and  reentry  phases  req>ectively 
against  the  missiles  from  a  given  launcher,  and  assume  that  exacdy  one  weapon  is  used  against 
each  in  each  phase.  If  X  survive  launch,  Wb  =  X  and  Wr  is  a  Binomial  random  variable  with  param¬ 
eters  X  and  q^.  Thus  EEW'r]  =  E[X],  and  E[Wr]  =  q^[X\. 

23  The  Final  Phase 

In  the  final  phase  the  probability  that  a  given  warhead  survives  an  attack  is  ^5.  Again  we  assume 
independence  among  aU  atten^ts  to  destroy  incoming  warheads.  Let  the  number  of  warheads  sur¬ 
viving  the  final  phase  from  the  i-th  incoming  missile  be  Z„  i  =  1, 2, ...,  Y.  Each  Z;  is  a  Binomial 
random  variable  with  parameters  n  and  <75,  so  F[Z,1  =  nqs  and  Var\Z'A  ~  die  number 

of  warheads  surviving  the  final  phase  (per  launcher)  be  //,  so 

H  =  iz, . 

11=  1 

Conditioning  on  7,  E[H\Y\  =  nYqs  and  Var[mY\  =  yVarEZJ  =  nYqs{\-qs).  Unconditioning, 


E[H\=nqiEm  (8) 

and 

Var[H\  =  «q5(l-95)£TIl  +  (9) 

where  £[7]  and  Var[7]  are  given  by  Equations  (5)  and  (6)  respectively. 

The  pmf  of  H,  pi^k),  is  found  in  a  similar  way  by  first  conditioning  on  7  If  7  =  0  (no  missiles 
survive  through  the  reentry  phase)  no  warheads  can  reach  the  target  area,  so  K))  =  1.  If 
7=y  >  0,  is  the  sum  of  j  identically  distributed  binomials  so  tiiatpaKk  \j)  =  b{k,  nj,  qs). 

Unconditioning, 


5 


(10) 


m 

J~k 

where  the  py0)’s  are  given  in  Equation  (7). 

Let  W^F  be  the  number  of  weapons  used  in  the  final  phase.  If  Y  missiles  survive  the  reentry  phase 
and  each  carries  n  warheads,  then  Wf  =  nY.  Thus  the  results  on  Y  can  be  used  to  calculate  the  mea¬ 
sure  of  interest  on  Wp 

Figure  2  demonstrates  the  model  by  showing  the  cumulative  tail  distribution  of  H  for  three  dif¬ 
ferent  sets  of  survival  probabilities.  For  aU  three  cases  the  number  of  missiles  per  launcher  (m)  is 
20,  and  the  number  of  warheads  per  missile  («)  is  10.  The  rightmost  curve  is  obtained  using  no  (or 
completely  ineffective)  counter  force  (qi  =  ^2  =  1)»  boost  and  reentry  survival  probabilities  (q^  and 
q^)  of  0.7,  and  a  final  phase  warhead  survival  probability  (qs)  of  0.4.  The  center  curve  is  obtained 
by  decreasing  q3  and  q^  from  0.7  to  0.6,  and  qs  from  0.4  to  0.3.  The  leftmost  curve  is  obtained  using 
the  original  set  of  parameters  but  decreasing  both  qi  and  qi  from  1  to  0.9.  Clearly  a  modest  increase 
in  kill  probability  in  counterforce  operations  from  0  to  0.1  has  a  dramatic  effect  on  the  number  of 
warheads  reaching  the  target  area.  An  increase  in  kill  probability  from  0  to  0.1  in  the  two  phases 
of  the  launcher  shows  a  drop  in  the  10-th  percentile  from  52  waiheads  to  23,  compared  to  a  drop 
from  53  to  3 1  for  a  similar  increase  in  kill  probability  in  the  boost,  reentry  and  final  phases.  Another 
way  to  interpret  the  three  curves  is  to  note  that  the  chance  of  at  most  20  waiheads  (10%  of  a  poten¬ 
tial  of  200)  reaching  the  target  area  is  3%  for  the  base  case.  With  a  given  improvement  in  active 
defense  this  increases  to  46%,  but  if  that  improvement  were  made  in  counterforce  instead  of  active 
defense  it  would  increase  to  87%.  These  numbers  are  shown  in  Column  2  of  Table  1.  Columns  3 
through  6  show  the  expected  number  of  weapons  used  in  each  phase.  A  small  improvement  in 
counterforce  effectiveness  sharply  decreases  the  number  of  weapons  required  for  active  defense. 
Note  that  the  zero  entries  in  columns  3  and  4  result  from  the  assumption  that  when  ^1=^2=  1»  it 
is  assumed  that  no  counterforce  is  attempted. 

The  next  section  contains  a  more  detailed  analysis  of  the  model  as  parameter  values  are  varied. 

3  Model  Analysis 

Throughout  this  section  results  are  demonstrated  using  kill  probabilities  ratho*  than  survival 
probabilities  p\  through  ps,  where  Pi  =  1  -  q,.  We  refer  to  a  kill  probability  vector  which  is  defined 
to  be  (pi,  P2.P3.P4.  Ps)-  For  example  (0, 02, 0.3, 0.5, 0.6)  represents  no  chance  of  killing  the 
launcher  in  its  outward  journey  to  die  launch  site,  a  20%  chance  of  kill  on  its  return  journey  to  re¬ 
load,  a  30%  chance  of  killing  the  missile  in  its  boost  phase,  a  50%  chance  in  its  reentry  phase,  and 
a  60%  chance  of  killing  each  warhead  in  the  final  phase. 

Theater  anti-missile  defense  today  consists  primarily  of  the  use  of  the  PATRIOT  system  in  the 
final  phase.  The  navy  Aegis  ship  anti-missile  ^fense  system  is  currently  being  considered  for  mod¬ 
ification  for  the  reentry  phase  of  anti-TBM  mission,  and  the  army  is  developing  the  THAAD  (the¬ 
ater  high  altitude  air  defense)  system  for  this  same  phase.  The  air  force  is  currently  developing 
boost  phase  systems.  Although  some  work  has  been  done  on  detecting  and  destroying  launchers 
prior  to  or  after  a  launch,  operational  experience  in  Desert  Storm  showed  that  current  systems  and 
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Figure  2:  Cumulative  Tail  Distributions  of  Warheads  Reaching  Target  Area 


Case 

Pr[H^} 

E[Wj^ 

Eiwy 

E[W^ 

E[Wr] 

9i  =  92=1.93  =  94®0.7,ft  =0.4 

0.03 

0 

0 

20 

14 

98 

9i  =  ?2=  1,  ft  *94  =  0.6,  ft  =  0.3 

0.46 

0 

0 

20 

12 

72 

9i  =  92  =  0.9, 93  =  94  =  0.7, 9,  =  0.4 

0.87 

5.2 

mm 

A.1 

3.3 

23 

Table  1:  Sample  Output  for  Numerical  Example 


operational  doctrine  are  ineffective.  This  current  state  can  be  modeled  by  setting  pj,  P2.  Pi 
all  equal  to  0.  We  can  set  Ps  at  some  value  depending  on  how  well  one  believes  the  PATRIOT 
works.  As  a  base  case  by  which  to  measure  possible  system  improvement  we  set  p^  to  0.7.  Thus 

Base  Case  Kill  Probability  Vector  =  (0, 0, 0, 0, 0.7).  (1 1) 
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Also  as  a  base  case  we  assume  that  a  launcher  can  launch  at  most  20  missiles  before  requiring  ma¬ 
jor  overhaul,  or  before  it  runs  out  of  missiles,  so  wi  =  20. 

We  look  at  three  measures  of  effectiveness  for  the  (random)  number  of  warheads  arriving  in  the 
target  area,  H.  These  are  (i)  the  mean  E\H\,  (ii)  the  median,  or  that  value  h  such  that  Pr{H^}  =  0.5, 
and  Oii)  the  ninetieth  percentile,  or  that  value  h  such  that  Pr[H^}  =  0.90.  We  also  look  at  the  ex¬ 
pected  number  of  active  defense  we£q>ons  required  (£[Wb],  fTWjJ,  and  fiTWp]),  and  the  expected 
number  of  count^force  weapons  (£tWc]).  We  first  look  at  today’s  case  where  there  is  only  one 
warhead  pa*  missile  (n  =  1),  and  show  how  some  performance  measures  are  affected  by  improving 
kill  probabilities  in  each  of  the  first  four  phases.  This  is  followed  by  a  similar  analysis  when  mul¬ 
tiple  warheads  are  considered. 

3.1  Single  Warhead  Analysis 

The  mean  numbers  of  warheads  (and  hence  missiles  since  we  are  assuming  one  warhead  per 
missile)  that  arrive  in  the  target  area  shown  plotted  in  Rgure  3  as  a  function  of  the  kill  probability 
at  a  particular  stage.  The  figure  contains  three  curves.  All  three  start  at  the  sante  point  (0,6)  because 
the  expected  number  of  warheads  reaching  the  target  area,  E[H],  is  6  when  m  =  20,  n  =  1,  the  base 
case  probabilities  are  given  in  (11),  and  Equations  (1),  (5),  and  (8)  are  used.  We  investigate  the  ef- 


Flgure  3:Mean  Number  of  Warheads  Reaching  Target 


feet  on  E[H\  of  increasing  each  of  the  four  zero  kill  probabilities  in  (11)  one  at  a  time. 
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The  upper  curve  is  found  by  increasing  the  kill  probability  of  either  the  boost  (p^)  or  reentry  (p^) 
phase  from  its  base  value  of  0  up  to  0.8.  In  either  case  it  decreases  linearly  with  a  slope  of  -6.  The 
middle  and  lower  curves  are  obtained  by  increasing p2  and Pi  respectively  ovct  the  same  range.  The 
difference  in  the  effect  of  a  small  increase  in  kill  probability  in  the  counteiforce  phases  when  com¬ 
pared  to  the  active  defense  stages  is  dramatic;  an  increase  from  0  to  0. 1  in  either  to  boost  or  reentry 
phases  reduces  £[//]  from  6  to  5.4,  whereas  this  same  increase  in  the  either  of  the  counteiforce  stag¬ 
es  reduces  it  fixim  6  to  approximately  2.5.  This  significant  improvement  is  caused  by  the  fact  that 
once  a  launcher  (and  its  crew)  is  destroyed  it  can  no  longer  fire  missiles,  causing  a  geometric  re¬ 
duction  in  E{H\.  In  the  active  defense  stages  a  kill  results  in  the  destruction  of  only  one  missile. 
The  small  improvement  in  increasing  p^  rather  than  pi  is  caused  by  the  fact  that  keeping  pi  at  zero 
means  the  first  missile  from  a  launcher  will  be  launched  for  certain,  whereas  increasing  p^  gives  a 
chance  to  destroy  the  launcher  before  its  first  missile  flies. 

Rgure  4  contains  a  similar  analysis  using  the  median  number  of  waibeads  reaching  the  target 
area  rather  than  the  mean.  Similar  results  are  found.  For  the  base  case  the  median  of  H  is  5.4.  In¬ 
creasing  the  boost  or  reentry  kill  probabilities  from  0  to  0.1  reduces  this  to  4.8,  whereas  this  in¬ 
crease  in  Pi  or  p2  reduces  it  to  1.3  and  1.6  respectively.  In  other  words,  using  a  kill  probability 
vector  (0.1, 0, 0, 0, 0.7)  there  is  a  fifty  percent  chance  that  fewer  than  1.3  warheads  will  reach  the 
target  area,  whereas  using  (0, 0, 0, 0.1, 0.7)  or  (0, 0, 0.1, 0, 0.7)  this  number  is  4.8. 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8 

Kill  Probability 

Figure  4:  Median  Number  of  Warheads  Reaching  Ihrget 

Figure  5  contains  a  similar  analysis  using  the  ninetieth  percentile  of  the  number  of  warheads 
reaching  the  target  For  the  base  case  there  is  a  ninety  percent  chance  that  the  number  of  warheads 
leaching  the  target  area  from  a  given  launcher  is  no  more  than  8.2.  Increasing  the  boost  or  reentry 


9 


Figure  5:  Ninetieth  Percentile  of  the  Number  of  Warheads  Reaching  Target 

kill  probabilities  from  0  to  0.1  reduces  this  to  7.6  whereas  an  increase  from  0  to  0. 1  in  P\Oipi  re¬ 
duces  it  to  5.5  or  5.7  respectively.  Although  by  using  this  measure  of  effectiveness  there  is  less  of 
a  difference  between  improving  counterforce  and  active  defense,  the  difference  is  still  significant 

We  now  turn  to  measuring  the  effects  of  changing  kill  probabilities  on  the  expected  numbers  of 
weapons  used  in  each  phase.  Starting  from  the  base  case  we  assume  that  a  zero  kill  probability  in 
a  given  phase  indicates  that  no  attempt  is  being  made  to  kill  the  launcher  or  missile  in  that  phase. 
Table  2  demonstrates  typical  results  that  can  be  obtained  fix)m  the  model  For  the  base  case  the 


Kill  Probability  Vector 

EmO 

E[Wj^ 

E[Wj,] 

E[Wj,] 

E\W^] 

Expected 

Warheads 

Killed/Weapon 

(0, 0, 0, 0, 0.7)-Base  Case 

0 

0 

0 

0 

20 

0.70 

(0,0,0,02,0.7) 

0 

0 

0 

20 

16 

0.42 

(0,0, 0.2, 0,0.7) 

0 

0 

20 

0 

16 

0.42 

(0,0.2, 0,0, 0.7) 

0 

4.94 

0 

4.94 

1.88 

(0.2, 0, 0, 0, 0.7) 

4.94 

0 

0 

0 

3.95 

2.11 

Table  2:  Effect  of  Increasing  Kill  Probabilities  on  Weapons  Numbers  and  Effectiveness 
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expected  number  of  weapons  used  per  launcher  when  no  attempt  is  made  to  destroy  the  missile  be¬ 
fore  the  final  phase,  and  assuming  one  weapon  for  each  warhead,  is  equal  to  the  number  of  missiles 
time  warheads  per  missile  that  a  launcher  can  launch.  In  this  example  that  is  20.  Also  for  flie  base 
case  the  expected  number  of  warhead  kills  per  weapon  is  equal  to  the  final  phase  kill  probability 
as  should  be  expected.  The  remaining  rows  in  Table  2  show  the  effect  of  increase  the  kill  probabil¬ 
ity  of  each  phase  in  turn  from  0  to  0.2.  Note  the  dramatic  drop  in  the  requirement  for  weapons  in 
the  final  phase  by  having  a  modest  effectiveness  in  counterforce  versus  the  same  effectiveness  in 
the  boost  or  reentry  phases,  hi  those  phases  a  modest  kill  probability  significantly  increases  the 
warhead  kills/weapons  used  ratio. 

3,2  Multiple  Warhead  Analysis 

We  repeat  the  analysis  of  Section  3.1  using  the  same  base  case  kill  probability  vector  shown  in 
(11)  and  twenty  missiles  per  launcher  (m  =  20),  but  in  this  section  we  assume  each  missile  carries 
ten  warheads  (n  =  10).  The  same  types  of  results  are  illustrated  in  Figures  6, 7  and  8  as  were  seen 
in  Figures  3, 4  and  5.  In  fact  since  the  mean  is  linear  in  n  the  curves  in  Figure  5  are  the  same  as 


Figure  6:  Mean  Number  of  Warheads  Reaching  Target,  Ten  Warheads  per  Missile 


those  in  Figure  3  except  the  vertical  scale  has  changed  by  a  factor  of  10.  There  is  no  simple  rela¬ 
tionship  between  the  median  or  the  ninetieth  percentile  and  n,  although  over  some  of  the  range  of 
the  kill  probability  the  relationship  appears  to  be  approximately  linear.  For  example,  from  Figure  4 
with  n  =  1  we  see  that  a  median  number  2  for  H  (90%  kill  of  the  twenty  possible  warheads)  can  be 
achieved  if  pi  or  P2  are  close  to  0.08,  whereas  in  the  boost  or  reentry  phases  we  would  pa  orp4  to 
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Figure  7;  Median  Number  of  Warheads  Reaching  Target,  Ten  Warheads  per  Missile 

be  0.56  to  achieve  this  success.  From  Figure  5  with  n  =  10  we  see  that  a  median  number  20  for  H 
(90%  kill  of  the  two  hundred  possible  warheads)  can  be  achieved  if  Pi  orpi  are  close  to  0.1,  where¬ 
as  in  the  boost  or  reentry  phases  we  would  P3  or  P4  to  be  0.65.  Similarly,  from  Figure  5  we  see  that 
to  achieve  a  ninetieth  percentile  of  2  when  n  =  l  requires  either  a  Pi  or  P2  of  about  0.28  or  a  ps  or 
P4  of  0.81;  from  Figure  5  a  ninetieth  percentile  of  20  when  n  =  10  requires  either  a  orpj  of  about 
0.30  or  a  pj  or  P4  of  0.79. 

Table  3  shows  the  expected  numbers  of  weapons  required  at  each  stage  and  the  expected  war¬ 
head  kills  per  weapon  when  n  =  10.  By  comparing  the  results  with  those  in  Table  2  it  is  clear  that 
the  required  expected  numbers  of  weapons  at  the  counterforce,  boost,  or  reentry  phases  does  not 
change  when  warheads  per  missile  increase  from  1  to  10,  but  the  number  of  weapons  in  the  final 
stage  increases  by  a  factor  of  ten.  These  results  should  be  expected  since  a  successful  kill  at  any 
phase  before  the  warheads  separate  is  assume  to  kill  all  n  warheads.  Note  that  the  expected  number 
of  warheads  irillml  per  weapon  increases  significantly  as  n  increases  the  earlier  one  can  attack  flie 
TBM  operation.  In  other  words,  counterfoice  is  increasingly  effective  as  the  number  of  warheads 
carried  by  the  missile  increases. 

33  Normal  Approximations 

For  given  values  of  m,  n,  and  a  kill  probability  vector,  it  is  easy  to  calculate  the  expected  value 
of  H  using  Equations  (1),  (5),  and  (8);  likewise  one  can  easily  find  the  variance  using  Equations 
(1),  (2),  (5),  (6),  and  (9).  But  to  find  percentiles  such  as  the  median  or  the  ninetieth  percentile  re- 
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Figure  8:  Ninetieth  Percentile  of  Warheads  Reaching  Target,  Ten  Warheads  per  Missile 


Kill  Probability  Vector 

E[Wm} 

E[Wj^ 

EWb) 

£[Wr] 

E[Wf\ 

Expected 

Wa^eads 

Killed/Weapon 

(0, 0, 0, 0, 0.7)-Base  Case 

0 

0 

0 

0 

200 

0.70 

(0, 0, 0,0.2, 0.7) 

0 

0 

0 

20 

160 

0.84 

(0, 0, 0.2, 0, 0.7) 

0 

0 

20 

0 

160 

0.84 

(0, 0.2, 0,0, 0.7) 

0 

4.94 

0 

0 

49.4 

3.41 

(0.2, 0, 0, 0, 0.7) 

4.94 

0 

0 

0 

39.5 

4.23 

Table  3:  Expected  Weapons  Numbers  and  Effectiveness  with  Ten  Warheads  per  Missile 


quires  the  distribution  function  of  H,  a  much  more  complex  calculation  using  Equations  (4),  (7), 
and  (10).  These  equations  were  used  to  find  the  curves  in  Figures  2, 4, 5, 7,  and  8.  Recall  that  H  is 
not  a  sinq>le  sum  of  independent  random  variables,  but  results  from  a  complex  set  of  five  random 
events,  the  first  two  of  which  have  a  truncated  geometric  distribution,  the  next  two  a  conditional 
binomial  distribution,  and  the  last  is  a  random  sum  of  these  weighted  binomials.  Even  so,  one 
might  suspect  that  its  distribution  is  approximately  normal  for  at  least  some  range  of  the  parameter 
values,  in  which  case  the  percentiles  can  be  estimated  using  only  the  mean  and  variance  of  H.  We 
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investigate  the  appropriateness  of  a  normal  approximation  for  the  median  and  ninetieth  po^ntiles 
of  if  in  this  section. 

Since  the  normal  is  a  symmetric  distribution  its  mean  and  median  are  equaL  Table  4  contains 


Kill  Probability  Vector 

Ten  Wariieads  per  Missile 
(n  =  10) 

Chie  Waihead  per  Missile 
(«  =  1) 

Median 

Normal 

Approximation 

Median 

Normal 

Approximation 

(0, 0, 0, 0, 0.7)-Base  Case 

59.3 

60.0 

5.4 

6.0 

(0,0, 0,0.2, 0.7) 
or 

(0,0, 0.2,0, 0.7) 

47.6 

48.0 

4.2 

4.8 

(0, 0.2, 0, 0, 0.7) 

10.3 

14.8 

0.6 

1.5 

(0.2, 0, 0, 0, 0.7) 

7.1 

11.9 

0.2 

1.2 

(0.2, 0.2, 0.3, 0.5, 0.7) 

NA 

2.3 

NA 

0.2 

Table  4:  Normal  Approximation  for  the  Median 


actual  medians  and  normal  approximations  for  the  base  case  and  kill  probability  vectors  used  in 
Sections  3.1  and  3.2,  and  an  example  that  assumes  positive  kUl  probabilities  in  aU  five  stages.  The 
normal  approximation  seems  to  perform  reasonably  well  for  the  ten  waihead  case  when  there  are 
zero  kill  probabilities  in  the  counterforce  stages;  it  does  less  well  in  the  single  warhead  case.  When 
Pi  and/or  P2  are  sigitificantly  larger  than  zero,  the  distribution  of  //  is  highly  skewed  and  the  normal 
approximation  for  the  median  is  poor.  The  entries  NA  (not  applicable)  in  the  table  indicate  that  the 
probability  that  H  is  zero  is  larger  than  0.5  so  that  no  median  value  exists. 

Figure  9  contains  cumulative  tail  distributions  (solid  lines)  and  normal  approximations  (dashed 
lines)  for  the  kill  probability  vectors  in  Table  3  and  one  warhead  per  missile.  For  none  of  the  ex¬ 
amples  is  the  normal  approximation  close  to  tiie  actual  distribution  excq)t  in  the  extreme  tails.  It 
is  particularly  poor  when  there  is  a  positive  probability  of  kill  by  counterforce. 

Figure  10  contains  cumulative  tail  distributions  (solid  lines)  and  normal  approximations 
(dashed  lines)  for  the  kill  probability  vectors  in  Table  3  and  ten  warheads  per  missile.  When  there 
is  no  counterforce  the  normal  approximation  is  close  to  the  actual  distribution  over  the  whole 
range,  but  again  there  are  significant  differences  when  there  is  a  positive  probability  of  kill  by  coun- 
terforce. 

As  one  might  expect  the  approximation  does  quite  well  when  H  is  &  fixed  (non-random)  sum  of 
binomial  random  variables.  Since  this  number  is  considerably  larger  when  multiple  warheads  are 
present  it  does  significantly  better  in  this  case.  With  positive  counterforce  probabilities  tiie  truncat- 
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Figure  9:  Cumulative  Tail  Distributions  and  Normal  Approximations,  n-1 

geometric  distribution  of  the  number  of  missiles  launch  leads  to  skewing  of  the  distribution  of 
In  this  case  the  normal  approximation  shows  significant  error. 

It  is  not  recommended  that  the  normal  approximation  be  used  for  the  median  (or  equivalently 
tiiat  the  median  and  mean  be  assumed  to  take  on  the  same  value).  Nor  is  it  recommended  that  it  be 
used  as  an  approximation  to  the  tail  distribution  unless  multiple  warheads  are  assumed  to  be 
present  and  the  only  significant  source  of  uncertainty  is  in  the  final  stages  of  the  TBM  operation. 


.5 


Warheads 

Figure  10:  Cumulative  Tail  Distributions  and  Normal  Approximations,  n  =  10 


4  Conclusions 

The  model  in  this  report  shoves  that  both  counteiforce  and  active  defense  will  form  essential 
parts  of  any  future  successful  system  for  theater  ballistic  missile  defense.  Without  counterforce  it 
will  be  relatively  easy  for  the  enemy  to  overwhelm  a  feasible  active  defense  system.  A  system  diat 
can  successfully  destroy  launchers  and  their  crews  will  provide  considerable  leverage  in  reducing 
die  numbers  of  active  defense  weapons  required;  this  leverage  increases  dramatically  as  the  num¬ 
ber  of  warheads  on  each  missile  increases.  The  model  allows  the  calculation  of  percentiles  of  the 
numbers  of  warheads  destroyed  rather  than  simple  expected  values. 

Past  eiqierience  in  finding  and  destroying  launchers  has  demonstrated  little  success  in  this  areas. 
As  was  discussed  in  Conner,  Ehlers,  and  Marshall  [1993],  success  will  most  likely  require  a  far 
more  structured  tqiproach  than  has  been  used.  A  model  for  such  a  structure  is  that  used  in  anti-sub- 
maiine  warfare  where  great  experience  has  been  gained  in  the  past  fifty  years  at  finding  and  de- 
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straying  toipedo  underwater  missile  launchers.  It  is  expected  that  the  successful  counterforce 
against  launchers  on  land  will  require  efforts  in  cuing,  search,  detection,  localization,  classification 
and  destruction.  Current  efforts  can  be  thought  of  as  attempting  to  skip  from  cuing  (for  example, 
flaming  datum  information  after  launch)  to  attack.  Future  reports  will  consider  how  one  might  best 
accomplish  the  in-between  phases  to  produce  successful  counterforce  against  mobile  missile 
launchers. 
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APPENDIX 


Expressions  are  derived  below  for  the  first  two  of  a  non-negative  integer-valued  random  variable 
in  terms  of  its  cumulative  tail  distribution.  Let  X  be  such  a  random  variable  and  letpi  =  Pr{X  =  i} , 
1  =  0, 1,2,...  . 

First  Moment 

£1X]  =Pi  +  ^2  3p3  +  4p4  + ... 


=P1 

+  P2 

+  P2 

+  P3 

+P3 

+  P3 

+  P4 

+  P4 

+  P4 

+P4 

+ 

=Pr{X>0 }  +  Pr{X>l }  +  Pr{X>2}  +  Pr{X>3 )  +  . . . 


Result  1:  E[X]  =  5;Pr{X>i} 

i-O 

Second  Moment 

£[X^]  -Pi  +  4p2  +  9p3  +  16p4  +  . . . 

=Pi 

+  P2  +P2  +P2  +P2 

+  P3  +  P3  +  P3+  P3  +  P3  +  P3+  P3  +  P3  +  Pi 

+  P4  +P4  +P4  +P4  +P4  +P4  +P4  +P4  +P4  +P4  +P4  +P4  +P4  +P4  +P4  +P4 
+  ... 

=Pr{X>0}  +  3Pr{X>l}  +  5Pr{X>2}  +7Pr{X>3}  + ... 

Result!:  £[X^]  =  (2i  + l)Pr{X>i}  =  2^JPr {X> i}  +£[X] 

i-O 


19 


INITIAL  DISTRIBUTION  LIST 


Defense  Technical  Infonnation  Center 
Cameron  Station 
Alexandria,  VA  22304-6145 

library,  Code  52 
Naval  Postgraduate  School 
Monterey,  CA  93943-5002 

Research  Office,  Code  08 
Naval  Postgraduate  School 
Monterey,  CA  93943-5002 

Department  of  Operations  Research,  Code  OR 
Naval  Postgraduate  School 
Monterey,  CA  93943-5002 

COL  R.  W.  Grayson,  USAF  Ret, 

POET 

1225  Jefferson  Davis  Highway,  Suite  300 
Arlington,  VA  22202 

Rich  Haver 

Executive  Director  for  Litelligence  Community  Affairs 
RM  5E56,  Original  Headquarters  Building 
Central  Intelligence  Agency  Headquarters 
Washington,  DC  20505 

Mr  Richard  Sokol 

Ballistic  Missile  Defense  Office 

Washington,  D.C.  20301 

David  Mosher 
National  Security  Division 
Congressional  Budget  Office 
Second  and  D  Streets,  S.W. 

Washington,  D.C.  20515 

NeilM,  Singer 
National  Security  Division 
Congressional  Budget  Office 
Second  and  D  Streets,  S.W. 

Washington,  D.C  20515 

R.  William  Thomas 
National  Security  Division 
Congressional  Budget  Office 
Second  and  D  Streets,  S.W. 

Washington,  D.C.  20M5 


Office  of  the  Chief  of  Naval  Operations 
N-85T  (I>r  Frank  Shoup) 

Washington,  D.C.  20350 

Office  of  the  Chief  of  Naval  Operations 
N-81 

Washington,  D.C  20350 

The  Joint  Staff 
Joint  Strategic  Planning  Staff 
Offutt  Air  Force  Base 
Nebraska,  68113-5001 

Professor  Michael  Bailey,  Code  OR/Ba 
Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Professor  James  N.  Eagle,  Code  OR/Er 
Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Professor  Donald  P,  Gaver,  Code  OR/Gv 
Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Professor  Wayne  P.  Hughes, Code  OR/HI 
Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Professor  Patricia  A.  Jacobs,  Code  OR/Jc 
Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Professor  Gordon  Nakagawa,  Code  OR/Na 
Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Professor  Samuel  H.  Parry,  Code  OR/Py 
Department  of  Operations  Research 
Naval  Postgraduate  School 
Montorey,  CA  93943 

Professor  Peter  Purdue,  Code  OR/Pd 
£>epartment  of  Opoations  Research 
Naval  Postgraduate  School 
Monterey,  CA  93943 


Professor  Michael  G.  Sovereign,  Code  OR/Sm 
Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Professor  Alan  R.  Washburn,  Code  OR/Ws 
Department  of  Opmtions  Research 
Naval  Postgraduate  School 
Monterey,  CA  93943 
Professor  Lyn  R.  Whitaker,  Code  OR/Wh 
Dq>artment  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  CA  93943 

CAPT  George  W.  ConnCT,  USN,  Code  OR/Co 
Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Professor  Kneale  T.  Marshall,  Code  OR/Mt 
Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  CA  93943 


